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Abstract 

We study hydrogen plasmas at mag netic fields B ~ 10^^ - lO" G, densities p ~ 10"^ - 10^ g cm and 
temperatures T ~ 10^"'' — 10^'^ K, typical of photospheres of middle-aged cooling neutron stars. We construct 
an analytical free energy model of the partially ionized plasma, including into consideration the decentred 
atomic states, which arise due to the thermal motion across the strong field. We show that these states, 
neglected in previous studies, may contribute appreciably into thermodynamics of the outer atmospheric 
layers at p < 1 g cm"'' and typical B and T. We take into account Coulomb non-ideality of the ionized 
component of the plasma affected by intense magnetic field. Ionization degree, occupancies and equation of 
state are calculated, and their dependences on the temperature, density and magnetic field are studied. 

1 Introduction 

Magnetic fields B ~ 10^^ — 10^"^ G typical of isolated neutron stars qualitatively modify many 
physical properties of matter ||l|, 0] . It was suggested that the outer layers of the neutron 
stars may be composed of hydrogen at temperatures T ~ 10^'^ — 10^'^ K Q. Thus the study 
of hydrogen plasmas at such B and T is of great practical importance for astrophysics. For 
studying the magnetized matter, Thomas-Fermi-like methods were used starting from 1970 
iQ] (see ref. Q for recent results and references). It is well known, however, that they are 
not well suited for light elements. Here we employ the free-energy minimization method. 

The motion of charged particles in a magnetic field is quantized into Landau orbitals. 
The magnetic field is called strongly quantizing if the free electrons populate mostly the 
ground Landau level |||]. This occurs when the electron cyclotron energy hujc = heB /{rriec) 
(where h, e, rUe and c are the Planck constant, electron charge, electron mass and speed of 
light, respectively) exceeds both the thermal energy ksT and the electron Fermi energy ep 
— that is for temperatures T <^Tb and densities p < pB, where 

h^B B 

Tb = 3.16 xW^^K, pB = 0.809 j^/' gcm~\ ^ ^ = (i) 

mice-^ 2.35 x 10^ G 

The atom in a strong magnetic field 7 ^ 1 is compressed in the transverse directions to 
the size of the "magnetic length": = {he/ eBY^"^ = ao7~^^^, where oq = h?/{mee^) 
is the Bohr radius. The ground-state binding energy grows logarithmically with B and 
exceeds the ground-state energy of the field-free atom by order of magnitude at i3 ~ 10^^ G 
||^]. Ionization equilibrium of atoms in strong magnetic fields has been first discussed in 
ref. However, that pioneering work neglected modifications of the atomic properties 
caused by the thermal motion of the atoms across the field. These motional modifications 
arise from the coupling between the centre-of-mass motion across the field and the relative 
electron-proton motion. These effects were appreciated by Ventura et al. 0, but quantum- 
mechanical calculations of binding energies and wave functions of hydrogen atoms in any 
states of motion in the strong magnetic fields have been carried out only recently |^. 
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Lai and Salpeter Q (see references therein for earher work) considered the ionization 
equihbrium of strongly magnetized hydrogen using a crude approximation for binding ener- 
gies of moving atoms which missed the so-called decentred states with a large electron-proton 
separation B. The same approximation was used in ref. |1C], devoted to the low-density 



equation of state. Here we employ new fitting formulae to atomic energies and sizes [11 
based on the previous study Q, valid for any state of atomic motion. We construct an 
analytic model of the plasma free energy and derive and solve a generalized Saha equation. 



2 Free Energy Model and Generalized Saha Equation 

We consider a plasma consisting of A'^e electrons, Np = protons, and Nh hydrogen 
atoms in a volume V, and write the Helmholtz free energy as -F = -Fid + Fex, where Fid = 

(e) (p) 

Fj^ + -|- is the sum of the ideal-gas free energies of the electrons, protons, and 
neutral species, respectively, and Fex is the excess free energy. 

For the ideal gas of electrons, the pressure and number density are 

where Ip{x) = /q°° t^dt/ (e*~^ -|- 1) is the Fermi integral, ^ IJ-e — NHlUc, fJ-e is the chemical 
potential, /? = (ksT)'^, Ag = h\/2-Kj3 /nie, gN>i = 2, and gN=o = 1- The free energy is 
given by F^^^'^ = ^e-^e — PeV, where is found using an algorithm described in ref. [p^]. 

In the strongly quantizing regime, the Fermi energy is ep = 27r^/i^ (o^ne)^/me, which 
differs from the non-magnetic case by a factor (4/3)2/3(p/p^)4/3. Thus the de generacy is 
strongly reduced at p <C PB- Furthermore, in the non-degenerate regime (fesT ^ e_p), we 



have F.f = NeksT [ln(27ra^A, 

For the protons, which are non-degenerate, we have 

/SFjf V^P = ln(27ra^Apnp) + In [1 - exp{-fihu;cp)] - 1, (3) 

where uJcp = {me/mp)uJc is the proton cyclotron frequency. Here, for sake of brevity, we 
drop the zero-point energy ^hwcp and the spin energy zizjgphucp, where gp = 5.585 is 
the proton spin gyromagnetic factor. These terms are the same for free and bound pro- 
tons. Taking them into account yields an additive contribution: AF = N(){hwcp/2 — 
kBT\n[2cosh{PgphuJcp/4:)]}, where No is the total number of protons (free and bound). 
AF does not affect ionization equilibrium and pressure. 

For the excess free energy of the ionized component, a general fitting formula in the 
non-magnetic case is given in ref. |]13[. It is known that thermodynamics of classical 
Coulomb plasmas is not affected by the magnetic field, which, however, affects the quantum- 
mechanical contributions to Fex. These effects have been studied only in the low-temperature 
or low-density regimes (e.g., ref. ||l^ and references therein). Here we use a scaling fr' 



eff 
s 



svs) of the density parameter = (47rneaQ/3) at a fixed Coulomb parameter F = 
j3e'^/{aQrs) in the formulae of ref. |13|. The scaling is devised so as to reproduce the low- 



density, high-temperature results presented in ref. ||10|, as well as other known limiting 
cases. For the contribution of electron-electron and electron-ion interactions in F^ex, the 
scaling factors are See = (1 + Om/Oo)/ [l + {9m/0o) exp(-6';;;^)/i] and Sie = 1//|, where 
^0 = 2 (9-7r/4)~^/^r5/r and 9m = 87^rf/(9-7r^r) are the non-magnetic and magnetic degen- 
eracy parameters, respectively, and the factors /i and /2 (depending on (3fiuJc) are given in 
ref. 0. 

The ideal-gas contribution of the magnetized atoms reads 

(^^if^ = Y.j ^^K^NsuiK^) {In [nH\]iwUK^)/Zj^ - l] , (4) 
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where s and v relate to electronic excitations, Ns,^ = {Xh /2ttH)'^ NhWsu^^'^'"' / are the 
atomic occupancies per unit phase space of the transverse component K± of the pseudo- 
momentum K which characterizes the atomic motion in the magnetic field, 'WsuiK±) and 
Xsiy{K±) are the occupation probabilities and binding energies of the moving atom, and 
Zw = {^H /'^'^Kf'Ylisy I ^ i-"^ su{K e'icgi[l3xsu{K _^_)\ is the internal partition function. 

The contribution of atoms in the nonideal part Fex of the free energy is calculated in 
the hard-sphere approximation using the van der Waals one-fluid model by analogy with 



ref. |14]. Its straightforward generalization to the magnetic case involves the composite 



quantum number k = {siyK±), so that includes now integration over K±. The hard- 



sphere diameters are set equal to the effective atomic sizes given in |11]. The occupation 
probabilities are then given by formulae derived in ref. extended to the magnetic case. 

Our model is valid as long as the formation of molecules may be neglected. In order to 
quantify the range of validity, we estimate the abundance of II2 molecules following ref. , 
but with inclusion of the non-ideal effects. 

Minimization of the free energy yields the ionization equilibrium (generalized Saha) 
equation: 

uh = npne(ApAe/A|^)(27ra^)^ [1 - exp{-(3huJcp)] Z^exp(A), (5) 

where A = f3fie — lii(27ra^Aene) + (3 dne/dlnrie — dPe/drif. takes into account effects of 
electron degeneracy and population of excited Landau levels. 



3 Results and Discussion 

Figure ^ shows selected results obtained for B = 10^^ G. The left panel shows the neutral 
fraction of atoms fn = Nh/Nq and molecules fH2 = '^Nh2/^o at T = 10^ K. For com- 
parison, we plot the fraction of atoms in the centred states, fn according to ref. and 
/h in the non-magnetic case. Long dashes display the fraction of atoms that satisfy the 
Inglis-Teller (IT) criterion and thus can be identified in optical spectra of the plasma. The 
IT fraction is estimated according to the formula ri^J- ~ exp[— np(4/K)^] (cf. Eq. (31) of 
ref. |l^). We can see that (a) the strong magnetic field increases the non-ionized fraction 
and shifts the region of pressure ionization to much higher p (compare the solid line and 
triangles in the left panel of Fig. |^), (b) the approximation of ref. reproduces only the 
abundance of the centred atoms at low density and fails at high density where the pressure- 
ionization effects are important, and (c) at the low density, the decentred atomic states are 
significantly populated. 

The right panel demonstrates the equation of state, which is seen to be much softer than 
(a) in the non-magnetic case (mainly because of the electron degeneracy "taken away" by 
the strongly quantizing field, but also due to the increased neutral fraction) and (b) in the 
magnetic but ideal proton-electron plasma (because the Coulomb interactions yield negative 
contribution to the pressure). 

The obtained results are used for modelling neutron-star atmospheres. In particular, 
the IT fraction of atoms, multiplied by the absorption cross sections calculated in ref. p5[| , 
determines an atomic contribution to atmospheric opacities. Preliminary calculations of 
the opacities, carried out with a simplified F^x, were presented in ref. The more 

elaborated model of the plasma described here confirms qualitative results of that work. An 
important conclusion is that the bound species contribute significantly to the absorption at 
B = 10^2 _ iQi3 g^gj^ relatively high T ~ 10^ K. 
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Figure 1: Left panel: Non-ionized fraction of atoms in any states (solid line), atoms in 
the centred states (short-dashed line), molecules (dash-dot line), and the weakly perturbed 
atoms contributing to the optics (long-dashed line), compared with the non-magnetic case 
(triangles) and the approximation |^ (dotted line). Right panel: Pressure isotherms of 
magnetized hydrogen plasma (solid lines) at lgr[K] = 5.5, 6.0 and 6.5, compared with the 
non-magnetic case (dashed lines) and with the ideal magnetized plasma (dotted lines). 
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